General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



FACILITY FORM 602 


X-66G-71-117 

PREPRINT 


i 


0.. ** ' T; 






IONIZATION AND HEATING 
OF THE GUM NEBULA 
BY ENERGETIC PARTICLES FROM 
THE VELA X SUPERNOVA 


R. RAMATY 
E. A. BOLDT 
S. A. COLGATE 
J. SILK 


’’A 


MARCH 1971 


t 



GODDARD SPACE FLIGHT CENTER 

RRFENBELT. MARYLAND 


N ^l - 2 1653 


^cce^o^w«mc'er) 

(THRU) 

IsJlIsEL 

(NASA CR OR TMX OR AD NUMBER) 

(CODE) 



(CATEGORY) 


IONIZATION AND HEATING OF THE GUM NEBULA BY ENERGETIC 
PARTICLES FROM THE VELA X SUPERNOVA 

by 

R. Ramaty and E. A. Boldt 
NASA/ Goddard Space Flight Center, Greenbelt, Md. 

S. A. Colgate 

New Mexico Technical Institute, S sco/ ro, N. M. 

J. Silk 

Astronomy Department, University of California, Berkeley 

ABSTRACT 

We investigate a model in which the Gum Nebula is ionized and 

heated by energetic particles from the supernova associated with the 

Vela X remnant. We find that, the observed ionization could be produced 
4 4 

in 10 to 2 x 10 years by the ejection of about 1 to 2 % of the pro- 
ducts of silicon burning at a mean energy of a few MeV/nucleon. We 
investigate the consequences of this model for the ionization and heat- 
ing of the interstellar medium, the generation of the light elements, 
X-ray production, and observable cosmic rays. 
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INTRODUCTION 

A physical model of the Gum Nebula has recently been presented by 

62 

Brandt et al. (1971) In which the nebula contains about 2 x 10 free 

electrons and is, assumed to be a fossil Stromgren sphere of the Vela X 

supernova. These authors indicate that if the ionization were produced 

by ultraviolet photons released by the supernova outburst at an average 

energy of 15 eV per ionizing photon, then the total energy in the ultra- 

51 

violet required to ionize the nebula is 5 x 10 ergs. In a more recent 
publication, Alexander et al. (1971) raise the possibility of a re* 
duction in size by up to a factor of 2, thereby reducing the total 
ionized matter by as much as a factor of 4 and reducing the total energy 

Cl 

of the ionizing photons to the range of 1 to 5 x 10 ergs. 

In this communication, we wish to discuss the alternative possibility, 

namely that the ionization may have been produced by fast charged particles 

released by the supernova. Since a fast particle of energy^ 0.1 MeV/nucleon 

expends some 36 eV to produce an ion pair in cold matter (Dalgarno and 

Griffing, 1958), the energetic particle output of Vela X must have been 

51 

at least 2 to 10 x 10 ergs in order to produce the ions in the Gum 
Nebula. This is not inconsistent with estimates of the total energy re- 
lease of supernovae. 

As we shall presently see, 

however, the required ionization and its spatial distribution depend on 
the assumed charge, mass and energy of the ionizing nuclei, on the age 
of the supernova aid on the matter distribution in the nebula, as well as 
on the propagation mode of the charged particles. 
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SUPERNOVA EJECTA 

The supernova phenomena is still only qualitatively understood and 
indeed | current controversy centers around the relative importance of 

i o 

the C thermonuclear detonation (Arnett 1969) versus neutron star pro- 
duction, (Colgate and White, 1966). In the present case, the formation 
of the neutron star cannot be disputed, but regardless, the ejected com- 
position ansi velocity distribution should not be drastically different 
for the two cases. The thermonuclear detonation of a carbon star (1.5 
i^) releases approximately 2 x 10 ergs. The formation of the Gum Nebula 
requires roughly 1 to 5 times this energy in ionization^ 
an amount which is certainly within the probable error of the 

explosion energy added by the dynamical formation of the neutron 
star. A fraction of the binding energy of the neutron star (*■ 10^ ergs) 
is presumably available to augment the thermonuclear detonation. 

The velocity distribution of the ejected matter is affected in the 
mean by the total explosive energy available, but the form is only slight- 
ly affected by the two processes. The simplest assumption is that the 

exploding gas expands as a uniform sphere (velocities calculated numerically are 
indeed in agreement with this, model). The presence of the neutron star perturbs 

the trajectories of a small inner mass fraction (Colgate 1971) and 
indeed the re-implosion of the central high density matter 
removes from observation what otherwise would be a prohibitive mass of 
neutron-rich matter. 

In addition to the uniform expansion approximation, we take note of 
the increasing strength of the explosion shock wave in the envelope of 
the star. Again numerical calculations show that this speed-up can be 
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roughly approximated by a power law In maus fraction, F, where F Is the 
fraction of the stellar mass external to the Lagrange coordinate posi- 
tion In question. 


Then: 
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The corresponding number spectrum of nuclei of mass number A in the non- 


relativlstic region is given by 

Me) • f ^ * f 


(£/£.f /2 ; e<e. 
(ec/^r ■> £>£< 


where m and Mej(A) are the proton mass and the total ejected mass in 
P 

nuclei of mass number A, respectively, and £ is kinetic energy per 
nucleon. The total energy output of the supernova in particles of mass 
number A depends on the transition energy £ 0 (which we assume to be the 
same for all nuclei) and is given by 


W = 4 £0 . 


Using this velocity distribution, Colgate and McKee (1969) showed that 
the light curve of a type I supernova could be created with a best fit 
to observations assuming 

1. The products of si lie on- burning resulted in the ejection of 
Ni^ and Si^ in a , ratio of 1 to 2 . 

2. An initial pre-supernova star ejects approximately 1/2 its 
mass composed of the above products. 


3. That the velocity v^ is chosen such that the heat given off by 
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Ni^— Co"***— Fe^ is re leafed by the progressive trans- 
parency of expansion resulting in the approximate conditions 
£©=“ (Mej/MD)^ MeV/nucleon and W 2 x 10^* (Maj^M^ ergs. 
Therefore, a small adjustment in the ejected mass in the range 1 to 2 
MZ) will cover the range of required ionization energies. The maximum 
I' possible energy is likely to be limited by the neutron star binding energy 
W * 1053 er g 8 or g 30 MeV/nucleon. In subsequent calculations we 
shall adopt £ 0 equal to 2, 10 and 20 MeV/nucleon. 


*■ THE NEBULA 

# 

' The time elapsed since the supernova explosion can be estimated 

* 

from the age of the pulsar PSR 0833-45, which is generally believed to 

W 

■-* 

t be the neutron star remnant of the Vela X supernova. This age can be 

W 

I determined from the present period and rate of increase of the period 

| and equals 11,500 years (Reichley, Downs and Morris, 1970). Shklovsky 

i 

•. (1970), however, suggests that the age of the supernova remnant as 

deduced from its radio surface brightness is about 30000 to 50000 years, 

A 

‘■the discrepancy can be understood in terms of repeated speedups of the 

* j 

Vela pulsar. He shall assume that t ■ 10 4 years is a reasonable lower 
%' limit (ft. the age of the nebula, with t < 5 x Hr years A possible upper 

I limit. 

% 

4 

According to Brandt et al. (1971) the comparison of the emission 
» measure of the Gum Nebula and the dispersion measure cf the pulsars in 
or behind the nebula Indicates that the matter is concentrated into 


and that 


clumps or clouds. Let (n e ) c and (n H ) c be the electron and hydrogen 
densities in the clouds and (n e ) c * (n e ) c + ^ n H^c* According to the 
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measurements summarised by Brandt et al. (1971), the present values of 

the spatial averages, < and <»l H ^are about 0.16 cm"\ 1.63 
-6 -3 

cm , and 0.062 cm respectively, so that <Vt o ^3£0*22. Under the 
assumption that all the matter is concentrated into the clouds, these 
averages can be written as 

where V is the nebular volume and V eff 18 the total volume of the clouds. 
The ratio X ■ V/V eff is then given by X ■ < ~ (oS • 

If the clouds are spherical objects with a random spatial dis- 
tribution, the mean distance l between clouds along an arbitrary dir- 
ection is & ^gjr > where r is the radius of a cloud. For a 
"standard cloud" (Spitzer, 1968) r * 7 pc so that^S£600 pc. The 
comparison of the emission and dispersion measures toward the Vela X 
pulsar indicates that the line of sight in that direction intersects 
at least one cloud. Since ( cannot be larger than the distance to the 
pulsar, f 5 pc. Obviously, both r and & may be highly variable 
throughout the nebula. 

The mode of propagation of the energetic particles in the Gum 

Nebula is essentially unknown. The high degree of isotropy ( 

11 

of the galactic cosmic rays with energies > 10 eV (Elliot, Thambyahpillai 
and Peacock, 1970) implies that the streaming velocity in interstellar 
space at these energies is less than ~ 20 km /sec or 2 x 10"^ pc/year. 

This velocity is obviously much too low to account for the ionization 
in < 5 x 10^ years of the entire GUm Nebula which appears to extend 
some 400 pc in the galactic disk (Brandt et al. 1971). The low anisotropy 
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of the galactic cotalc rayt hat recently been Interpreted by Llngtnf altar, 
Raaaty, and Flak (1971) in tarmt of particle propagation by compound 
dlffutlon. Such dlffuaion comblntt tht effecta of random walk of tht 
field llnea (Joklpli and Parker, 1969) with one-dimensional dlffuaion 
along field lines. The net spatial displacement of par tic let under- 
going compound dlffuaion it proportional to t^ at opposed to the depend- 

1/2 

ence proportional to t ' at encountered in ordinary dlffutlon. This 
leads to a low cosmic ray propagation velocity and can account for the 
isotropy observed at high energies. 

Since diffusive propagation along field lines it the result of 
scattering of the particles by magnetic irregularities, it is possible 
that in the few MeV/nucleon range the gyroradii are much smaller than 
the scale sire of the irregularities, hence particles will no longer 
be scattered and diffusive propagation will reduce to regular motion. 

In addition, shortly after the supernova explosion, the energy density 
in low energy particles is much larger than that in field irregularities, 
so that even if scatterings were important, the irregularities are ex- 
pected to be swept along the field lines by the particles. In both 
cases, the particles propagate along the field lines with a mean 
velocity c^ , where |& * 0.15 for £ * 10 MeV/nucleon. If the character- 
istic length of the random field distribution in the Gum Nebula is of 
the order of the mean distance between clouds, particles could propagate 
to the edge of the nebula with essentially rectilinear motion and be 
stopped or deflected only by collisions with clouds. With * 0.15, 
particles will traverse 400 pc in 8200 years, a time which is less than 
the lower limit on the age of the nebula. 
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Whereas tht motion of tho particles in the intorc loud medlia la 
essentially ragular motion Along tht fiald lint a, panatration into tha 
cloud* ia probably govarnad by wave -particle intaractiona. Aa tha 
particlaa atraam down a danaity gradiant thay product hydromagnatic 
wavaa which acattar tha particlaa and limit tha atraaming (Wantsal, 1969; 
Kularud, and Pearce, 1969). Tha limiting atraaming valoclty Ufa can be 
obtained by balancing tha growth rata and damping rata of tha wavaa, and 
la approximately given by (Kularud and Cesarsky 1971) 


w « ■ [' 1 * * Bfc ^ 7 >) ] 


Her<>r) J (5) 

Here la the Alfvan apeed; B ia the magnetic field in clouda; R(«jj) c 

la the damping rate which ia temperature dependent and equala 1.5 x 10*9 

(n H ) c aec”* at 100 K and about 1.5 x 10’® (njj) c aec“* at 5 x 10^ K; p 

ia the momentum at which the streaming ia evaluated; and n cr (>p) ia the 

denaity of energetic particle* with momenta greater than p. 

Consider the Initial interaction of the energetic particles with HI clouds 4 a* a. 

lower limit for the particle density in clouds, we take n ■ N 0 /V 

cr 

where N 0 ia the total particle output of the supernova and V is the 

total volume of the nebula.For N Q corresponding Si 3 ® + 0.25 MD Fe$® 

and V - 1.35 x 10 63 cm 3 (Brandt et al.,1971) we get % l + 

(j*fi] for neutral clouda at 100 K,and 

for ionized clouda at 5 x 10^ K. Let us assume that the Initially neutral 

clumps in the Gum Nebula are similar to ordinary HI clouds for which 

(r>-e) c - 0*03 (Hje liming, Gordon, and Gordon. 1969). For these 

clouds in the Gum Nebula (n H ) c 14 cm" 3 , so that for B » 3 ~ 42 km/ sec. 

For ionized clouds in the Gum Nebula (n ) ** 10 cm" 3 and (n M V 4 cm" 3 , no 

e c He 

that for 1 3 * 3 /tig, If— 226 km/ sec . If we assume that the average transit 


v? j 1 *£$£ I '•■ * ^ r 't* ' ! 
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dlstance across s cloud Is of the order 5 pc, the transit times in 
neutral and ionized clouds become about 10^ and 7 x 10 3 years, re- 
spectively. 

These times can be compared with the energy loss times in the 
clouds. The energy loss per nucleon of a nonrelativistic particle of 
charge Ze in partially ionized hydrogen is given by (Hayakawa and Kitao 
1956) 

^ = 1,H6 y to" '**jC f ( ^ 1 t,fr e ) c ] £ * *’ MeV/w.l. ac*' (6> 

where £ is in MeV/nucleon. The time in which a silicon nucleus of 10 MeV/ 
nucleon will lose half its energy in neutral or ionized clouds, is about 
2000 years and 640 years, respectively. Both these times" are much shorter 
than the brans! t times across clouds, so that once a particle penetrates 
a cloud, it will stop and deposite its energy on a time scale less than 
about 2000 years. 

We now consider the energy deposition and ionization. The total 
energy loss to both ionization and heating, E t , of particles with 
initial spectrum N 0 ( £' ) in a fixed time t is given by 

£* = hl'M A O' - 

% J c (7) 

where £(t,t) is the energy per nucleon as a function of time of a 
particle of initial energy per nucleon . Because both (n e ) c and (n H ) c 
are time dependent as a result of the varying degree of ionization of 
the clouds as a function of time, the solution of equation (6) has to be 


written as 

e(e‘,t) = £'[/- 
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where t If In years, £ and in MeV/ nucleon, and 

nvf+itne * £ [ <lt\ 

The energy loss to ionization alone is given by 

r _ ( U \c< W /£•)(*_(-."■’* — 4i Ji' 

^ * \ No '°i (n,) c ^(n t ) c dt 

where <Je/J t io given by equation (6). By introducing the average 
quantity 

, t 

/ — ^ — \ r (Y-e(*'.t)f f -- folk - 4t’ 

UnMrVej L ^ J ) 4t ’ 


(9) 


( 10 ) 


^ v ic* 

Che energy lose to ionlz stlon can be written as 


( 11 ) 


Ei * 


l ru 

\ n M +Ai 


•\ k 


htAi fie) ~ • (12) 

We have evaluated equation (12) with Ej. given by equation (7), N 0 by 
equation (2) and £ by equation (S). The results, for various values 
of £ a , are shown in Figure 1 as functions of £(17# + . The 

solid line is the energy loss of 0.5 MU Si^® 4- 0.25 MD Fe^® and the dashed 
line is the energy loss of 0.75 M3 H*. As can be seen, all curves as sump - 
totically approach the total energy W given by equation (3), 

6 2 

The total energy loss to ionization required to produce the 2 x 10 
ion pairs in the Gum Nebula is about 10®^ ergs. We neglect recombination 


in clouds, since for an electron temperature greater than about 5 x 10^ K 
(Alexander et al. 1971) the recombination time ( ~ 

Bates and Dalgarno, 1962) is greater than about 2 x 10^ years. If the 
particles lose their energy predominantly in neutral hydrogen 
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and the Ionization is produced by silicon and iron nuclei with So 
between 10 and 20 MeV/nucleon, tn„ must be about 4.5 x 10^ cm"** 
years. For * 14 cra“* , t * 3200 years. This is consistent with 

the total time spent by the particles in clouds. 

If the ionization is produced by protons with a spectrum of similar 
So an< * cm” ^ , t is about 24,000 years which exceeds the lower 

limit on the age as determined from the pulsar observations. Furthermore, 
as pointed out by Alexander et al. (1971) since only 1/65 of space in the 
nebula is occupied by matter, pressure equilibrium has not yet been es- 
tablished betweenthe clouds and the intercloud medium so that the ioniza- 
tion had to occur over a relatively short time. It is therefore more 
likely that the ionization and heating in the Gum Nebula have been produced 
by low energy silicon and iron nuclei on a time scale in clouds of less 

than 2000 years and not by low energy protons on a time scale of more 
4 

than 10 years . 

COMPARISON WITH RELATED OBSERVATIONS 

We now consider some of the consequences of this model regarding 
the heating and ionization of the interstellar medium, the generation 
of light elements, bservable X-ray production and the detectability of 
charged particles at earth. 

Several recent theoretical descriptions (Pikel'ner 1967; Balasubrah- 
manyan et al. 1968; Spitzer and Tomasko 1968; Spitzer and Scott 1969; 
Goldsmith, Habing, and Field 1969) of ionization and thermal equilibria 
for interstellar HI regions, based upon current models for the inter- 
stellar gas, indicate that the rate of ionization per hydrogen atom is 
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=- !0" 15 (sec H atom)" 1 . HJe liming, Gordon, end Gordon (1969) find 
that observed pulsar dispersion measures may be best fitted with such a 
model for ■ (2,5 + 0,5) x 10" 15 (sec H atom) -1 . For an equilibrium 
situation, the rate of electron-ion recombination is the most direct 
measure of ^ ; recent observations of H$ hydrogen line emission from 
interstellar HI (Reynolds 1971) yield a direct measure of the recombin- 
ation rate, and the corresponding ionization rate in the regions examined 

- 14 -i 

could be as high as 10 (sec H atom) . 

Since essentially all the gas in the Gum Nebula has been ionized by 
the supernova explosion, the average rate of ionization in interstellar 
hydrogen by such events would be given by 

< "5 > - 4sn (13) 

where V gN is the volume of the HII region formed by the supernova explo- 
sion and fg^ is the average supernova frequency per unit volume in the 
galaxy. If supernovae occur at a rate of 1 per 100 years in a volume of 
4 x 10 66 cm 3 , f gN =* 8 x l(f 77 cm‘ 3 sec" 1 . For V gN <. n e > ~ 2 x 10 62 elec^Or^ 
(Brandt et al. 1971)^we obtain n^^ 1*6 X 10" 14 sec" 1 cm"'*. Radio 

observations at 21 cm indicate a mean value for the neutral hydrogen 
density in the galactic plane of about 0,7 cm"’* (Spitzer 1968) j however, 
optical depth effects probably raise ^n^ to ~ 1 crn”^. A further in- 
crease by about a factor of 2 in the value of ^n^ may be due to the 
presence of molecular hydrogen in dense clouds. Hence the mean galactic 
value of as indicated by the parameters deduced for the Vela X super- 
nova is about 10" 14 sec" 1 (H atom)" 1 . This ionization rate could be 
reduced if the frequency of supernova explosions of the type that produced 
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the Gum Nebula Is lover than the assumed rate of 1 per 100 years or If 
the actual free electron content of the Gum Nebula is less than that 
estimated by Brandt et al. (1971). Furthermore, even if the average 
supernova output is similar to that of Vela X but occurs in a region 
with larger magnetic fields, so that the particles propagate over a 
snorter distance before they lose their energy, the ionized volume pro- 
duced by the supernova would be smaller than that of the Gum Nebula. 

In this case, the total ionization per supernova would be reduced since 
a larger fraction of the total energy would go into heating. Also, 
those supernovae which occur at more than about 100 pc away from the 
galactic plane (e.g. the Crab) will produce energetic particles which 
may escape from the disk without producing significant ionization. 

In order for HII regions such as the Gum Nebula to merge into the 
HI of the interstellar medium, the cooling time of the clouds should be 
much less than the time between supernova explosions in the volume of the 
nebula. For f gN * 8 x 10’ 77 cm“ 3 sec _1 and V gN * 1.35 x 10 63 cm 3 , this time 
equals 3 x 10^ years. However, as we shall see from the calculations of 
Cox and Tucker (1969), the cooling time in clouds could be several orders 
of magnitude shorter. For a gas of cosmic abundances in ionization equi- 
librium at 5 x 10^ K, the radiative loss rate is P =“ 2 x 10’ ^ 2 (n e ) c 
( n H) c ergs cnf^sec . The corresponding cooling time is t c =* 3/2 (n e ) c 
kT/P =“ 1700/ (n H ) c years. For (n H ) c = 4 cm" 3 , t c =* 420 years. Although 
this cooling time is much shorter than the age of the nebula, contemporary 
heating by cosmic rays could in fact maintain the observed high temperature 
(Alexander et al. 1971). In contrast, the rapid cooling would be in con- 
flict with ionization by ultraviolet photons which deposit their energy 
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in the first thousand years after the supernova explosion. However, if 
the abundances of the principal ions responsible for the cooling are 
depleted below their values corresponding to ionization equilibrium with 
n^j ■ 4 cm , the cooling time could be proportionally increased.^ 11 addition, 
if the dimensions of the clouds are sufficiently large, there is the 
possibility of the entrapment of the resonant radiation which leads to 
the cooling. This phenomenon may further increase the cooling time ob- 
tained from Cox and Tucker's (1969) calculation. 

We now consider the production of the light elements. According to 
Fowler, Reeves and Silk (1970) the average ionization rate of the inter- 
stellar gas by cosmic-ray protons of energies greater than 5 MeV/nucleon 
has to be less than about 10” (sec H atom)"*. A higher ionization rate 
by such cosmic rays would result in the production of lithium, beryllium, 
and boron, with abundances in excess of the observed values. However, 
if the interstellar gas is ionized by low energy heavy nuclei which have 
lower spallation cross sections for the production of Li, Be, and B than 
do protons and the medium nuclei (C,N.O), a much higher ionization i tie 
is permissible. 

28 

Consider the production of lithium by Si nuclei. The total number 
of Li nuclei produced by the Vela X supernova is o (n c ) c < t N 0 (>£^), 

where a is the production cross section, t is the time spent in clouds by 
the silicon nuclei with energies greater than the threshold energy £ ^ , 
and N q is given by equation (2) . Neglecting escape or destruction, the 
total lithium content of the galaxy equals the total production over 
its lifetime Tg, i.e., 

H u - <rK) c < c p (*/T) li (>e th ) r j 


( 14 ) 
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where T Is the average time between supernova explosions In the galaxy. 

The lower limit on the threshold for lithium production from Si^® on 
hydrogen is about 20 MeV/nucleon and the cross section at 150 Me V/ nucleon 
is about 2.5 mb (Shapiro and Silverberg 1970). From equation (2) with 
So “ 10 MeV/nucleon and * 20 MeV/nucleon we get N 0 (>£^) =“0.1 Mej/(Amp), 


ao 


For T ■ 10 AW years, <$> * 0.2 and Mej ■ 0.75 M2) we get that N^ =- 2 x 10*^(t/T). 

O 
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The hydrogen content of the interstellar medium is about 4 x 10 so that 
the ratio Li/H is 5 x 10“** (t/T). According to the data summarized by 
Reeves, Fowler, and Hoyle (1970), Li/H is less than about 10“^ so that t/T 
is less than 20. Since the time spent in clouds by silicon nuclei above 
threshold energy is on the order of the energy loss time, which for a 20 MeV/ 
nucleon particle may vary from about 1600 years to 5000 years, the lower 
limit on T ranges from about 80 years to 250 years. This result can be 
compared with the supernova rate required to produce the iron content of 
the galaxy. According to Clayton, Colgate and Fishman (1969), the galaxy 
contains about 3 x 10' MD ¥e J . If each supernova produces 0.25 M2) Ni J , 
which decay into iron, for Tg * 10^ years, T is about 80 years. We con- 
clude that the rate of production of light elements from supernova ejecta 
of the type that produced the Gum Nebula is not in conflict with the ob- 
served upper limits on the abundances of light elements. This conclusion 
is independent of whether the particles lose their energy in a predomin- 
antly neutral or ionized medium. 

X-ray emission is a necessary consequence of ionization and heating 
by cosmic rays (Boldt and Serlemitsos 1969). In clouds, the volume 
emissivity of bremsstrahlung X-rays arising from the collisions of ener- 
getic nuclei with ambient electrons (free or atomic) is given by 
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2 G© ^o) c cptri CT (£)A£. 

where cr 0 ■ 6.7 x 10’ 23 cm 2 , <y ■ 1/137, m e is the electron mass and n 


(15) 


Is the differential number density of the energetic particles in clouds. 

In order to estimate an upper limit to the X-ray production, we shall 
assume that all the energetic particles ejected by Vela X are currently 
In clouds at the characteristic energy £ 0 • The resulting X-ray energy 
spectrum is then essentially flat with an end point ho * (m&/mp) & a $ 
for m 20 MeV/nucleon, hv * 10 keV, 

According to the geometry of the Gum Nebula described by Brandt et 
al. (1971), we expect it to appear as a disk source of X-rays, extended 
in galactic longitude (aI =“ 90°) and relatively thin in galactic latitude 
(A# - 7°). Collecting radiation from an interval of length L penetrating 


the nebula, the intensity viewed by a detector of aperture (A 


i)< A £ 


would be expressed by 


il =■ nL 

AO nr L 
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Here ? * ?c Vgff/V is the average volume emissivity of the nebula and, 
with given by equation (15), may be expressed as 

i ' % > 

where <c|J> is the average velocity over the particle spectrum. For 
Mej = 0.5 MS Si 28 + 0.25 MS Fe 56 and « 0.2, * 1.2 x 10‘ 27 ergs 

cm' 3 sec’ 1 . Using a detector of (A^) 0 * 4°, viewing at $ (260 + 40)°, 

which spans the Gum Nebula, Cooke, Griffiths and Pounds (1969) observed 
an enhanced X-ray emission associated with the Hi 
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3 x 10 ergs (cm sec rad) above 1,4 keV, with an essentially flat 
energy spectrum extending to about 10 keV. The upper limits to this 
flux at energies above 12.5 keV, as set by Hudson, Peterson, and Schwart 2 
(1971), provide evidence that the spectrum Is characterized by a break 
in the vicinity of about 10 keV. Using the volume emisslvity as obtained 
from equation (17), L ■ 400 pc and (A ^) 0 * 4°, we get dl/d^ ■ 8 x 10*9 
ergs sec” A cm"*rad" . The fact that this upper limit is a few times 
larger than the observed X-ray flux means that a significant fraction of 
the particles ejected by the supernova have already lost their energy 
or that the characteristic energy £ e Is smaller than the value of 20 MeV/ 
nucleon that we have used in this calculation. 

In addition to bremsstrahlung, X-rays in the Gum Nebula should also 
be produced as line emission. These lines are produced by charge ex- 
change between the energetic nuclei and neutral H atoms, resulting in 
the capture of electrons to excited states. Cascades down to the ground 
state occur, which produce the analogue of L^c( and emission for the 

case of capture of K-electrons. Silk and Steigman (1969) have given a 
discussion of this process, with application to low energy galactic cosmic 


i f 


rays, 


If we assume that each nucleus captures two electrons in the life- 


time of the nebula and hence produces two X-ray photons, the time averaged 
volume emisslvity is 


/, v „ Me/ V>U 
= ^ Awy Vt 


( 18 ) 
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where V end t ere the volume end ege of the nebule, respectively, end 
hu ** 2 keV for silicon end hu “* 7 keV for Iron. For Mej * 0.5 MS) Si 3 ® 

+ 0.25 MD Fe 56 , t - 10 4 yeers end V - 1.63 x 10 63 cm 3 , I £ (2 keV) 

**■ 3.4 x 10*^® ergs cm” 3 sec - ^ and (7 keV) =“ 6 x 10"^® ergs cm”"*sec~^. 
Using equetlon (16), we find that the total X-ray flux in line emission 
is comparable to or somewhat greater then the total observed X-ray flux 
from the direction of the Gum Nebula. However , since charge exchange 
occurs only in neutral hydrogen, the contemporary X-ray flux in line 
emission could be smeller than the time-averaged flux that we have 

estimated. Since the energy of the nuclei at capture is 

not expected to be more than about 20 MeV/nucleon, the line broadening 
should be less than 20%. Therefore, detectable X-ray line emission 
could be observed with detectors of moderate energy resolution. 

Finally, we consider the implications of our model on direct cosmic 
ray observations at earth. The observability of low energy (£ 100 MeV/ 
nucleon) heavy nuclei from the supernova depends on whether these parti- 
cles escape from the nebula before stopping in a cloud and on a rather 
large but unknown amount of solar modulation which may include a sub- 
stantial energy loss in the interplanetary medium (Goldstein, Fisk, and 
Ramaty 1970). Above 100 MeV/nucleon, however, these nuclei are not 
expected to stop in clouds in a time which is short compared to the 
average time between supernova explosions. Therefore, it would appear 
that our model predicts silicon and iron intensities in the galactic 
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cosmlc rays which exceed their observed values (e.g. Meyer 1969). 

However, current Ideas concerning pre- supernova structure (Arnett 1970, 

S.A. Colgate, and G. LeFebre, In preparation) are compatible with a 
change in composition for mass fractions F < 10~^» For these mass 
fractions the composition could be consistent with that of the observed 
galactic cosmic rays. The energy corresponding to F * 10”^ is about 
500 MeV/nucleon for £ 0 m 10 MeV/nucleon. However, adiabatic energy 
losses associated with the expanding magnetic field in the vicinity 

of the supernova may decelerate these nuclei below ~ 100 MeV/nucleon, { 

•4 

where energy losses in the interstellar and/or interplanetary medium i 

prevent us from observing the particles. 


/ 
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SUMMARY 

The ionisation and heating of the Gum Nebula could be produced by 
the dissipation of the energy of about 1 to 2 M3 of the products of 
silicon burning ejected by the supernova Vela X on a time scale of the 

A 

order 10 years. The characteristic energy at injection of these parti- 
cles is on the order of a few Me V/ nucleon so that the total supernova 
energy in charged particles is a few times 10^ ergs. The propagation 
of these particles in the Gum Nebula is essentially regular motion along 
field lines so that the particles fill up the nebula on a time scale 
comparable to its age. For the complete ionisation of the nebula, it 
is essential that a majority of the particles should have stopped in 
clouds. 

The average rate of ionisation per hydrogen atom in the interstellar 
medium due to the formation of HIX regions similar to the Gum Nebula, at 
a rate of 1 per 100 years in the galaxy, is about 10” ^ (sec H atom)”'*'. 

This is independent of the nature of the ionizing radiation, (l.e., photons 
or charged particles). However, ionization and heating by charged particles 
allows a cooling time in clouds which is shorter than the age of the Gum 
Nebula. Contemporary heating by particles would result in X-ray production. 
The upper limit to the expected X-ray flux from the direction of the Gum 
Nebula exceeds the measurements by about a factor of 2 to 3. Therefore, 
even though there may still be continuous heating of the clouds in the 
nebula, a significant fraction of the initial energy has already been 
dissipated. 
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The amount of light nuclei produced during the Ionisation and heat- 
ing of the Gum Nebula Is consistent with the observed abundances of the 
light elements If supernovas similar to Vela X occur at a rate of about 
1 per 100 years In the galaxy. The heating of the Interstellar gas by 
low energy silicon and Iron nuclei is therefore not in conflict with the 
abundances of lithium, beryllium, and boron in the interstellar medium. 
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FIGURE CAPTION 


1. The total energy loss, E t , and the energy loss to ionization, E^, 
as a function of ♦‘he age and average gas densities in clouds. The 
quantities njj + 4 n e and (n^ + 4 n e )/n^ are defined in equations (9) 
and (11), respectively; £ o is the characteristic transition energy 
defined in equations (2) aid (3). 


